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QUASIL INEAR  SCATTERING  FROM  WAVES  DRIVEN  BY 
BEAM-PLASMA  INSTABILITIES 


INTRODUCTION 

Test  data  taken  during  the  Starfish  event  of  1962  displays  a 
marked  second  brightening  of  the  atmosphere  in  the  southern  conjugate 
region.  This  brightening  followed  the  initial  luminescence  in  the 
southern  conjugate  region  caused  by  the  deposition  of  the  initial  debris 
by  approximately  10  seconds. 

The  explanation  that  appears  most  probable  involves  scattering 
a  fraction  of  the  streaming  debris  out  of  the  loss  cone  as  the  debris 
proceeds  toward  the  southern  conjugate  region.  The  scattering  sources 
are  waves  excited  to  high  levels  driven  by  the  initial  fast  debris. 

The  scattered  debris  would  then  mirror  in  the  southern  conjugate  and 
stream  toward  the  magnetic  bubble.  Some  fraction  would  be  scattered 
back  into  the  loss  cone  region  by  the  loss  cone  instability.  The  remain¬ 
der  would  mirror  off  of  the  bubble  and  return  toward  the  southern 
conjugate  region.  During  the  return  the  loss  cone  region  could  be  filled 
by  loss  cone  instabilities  and  a  second  precipitation  of  debris  particles 
could  occur.  The  second  precipitation  would  involve  a  greater  proportion 
of  the  distribution  because  of  nonlinear  flattening  of  the  distribution 
function. 
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In  this  report  the  possibility  of  exciting  the  waves  which  are 

predicted  to  perform  the  scattering  is  addressed.  It  will  be  shown  by 

use  of  a  numerical  code  and  quasilinear  estimates  that  indeed  waves  can 

be  excited  and  will  grow  to  sufficient  levels  to  cause  quasilinear 

diffusion  of  the  streaming  debris  distribution  function. 

The  code  employed  is  a  linear  code  which  has  been  used  to  study 

beam-plasma  instabilities  in  tokamaks.  The  circumstances  of  the  problem 

are  similar  to  a  beam-plasma  interaction,  the  general  type  of  which,  have 

1-5 

been  studied  for  a  number  of  years  in  the  fusion  community.  The 

linear  code  calculates  the  energy  transfer  between  the  beam  and  waves  of 
interest  and  the  waves  with  the  background  plasma  to  determine  the  growth 
rates  and  energy  transferred  to  the  wave.  From  these  calculations 
quasilinear  estimates  are  made  to  determine  the  amount  of  scattering 
taking  place. 

There  are  a  large  number  of  waves  which  may  be  excited  by  a  beam 

of  high  energy  ions.  Among  them  are  the  ion  cyclotron  wave,  ion 

Bernstein  wave,  ion  acoustic  wave,  low-hybrid  waves  and  shear  Alfven 

waves.  Several  of  these  waves  are  of  high  frequency  and  candidates  for 

study.  However,  the  ion  cyclotron  wave  and  its  higher  harmonics,  usually 

called  the  ion  Bernstein  wave,  are  chosen  for  study  because  these  waves 

1  6-10 

have  been  predicted  theoretically  and  shown  experimentally  to  be 

excited  during  beam  plasma  interactions.  The  code  employed  for  this  work 

was  successful  in  predicting  quite  accurately  these  instabilities  in  a 

variety  of  plasma  devices  from  Q  machines,  to  Tokamaks  such  as  the  French 
10 

tokamak,  TFR,  and  PLT,  the  Princeton  Large  Torus. 
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In  the  following  section  the  details  for  calculating  the  linear 
growth  rates  for  the  waves  will  be  briefly  discussed.  The  quasilinear 


estimates  are  presented  in  the  following  section  and  finally  the 
conclusion  of  these  calculations  are  presented. 
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LINEAR  CALCULATION 


The  linear  calculation  is  based  on  a  perturbation  method  where 

4  11 

the  real  frequency  is  assumed  larger  than  the  imaginary,  0)^  >  ou^.  ’ 

The  code  calculates  the  rate  of  change  of  the  wave  energy  by  computing  a 
time  averaged  J_  •  E  for  each  species  given  by  the  following  formula 


=  <J  •  E> 


where  the  sum  j  is  over  the  plasma  species.  Physically  the  process  of 
energy  transfer  takes  place  through  resonant  processes  either  Landau 
cyclotron  or  anomalous  cyclotron  resonances.  The  beam  transfers  energy 
to  the  wave  and  background  plasma.  The  wave  in  turn  transfers  energy 
to  the  background  electron  or  ions.  If  the  beam  pumping  rate  is  larger 
than  the  damping  rates  the  wave  grows. 

Typically,  the  instability  is  parameterized  by  the  amount  of 
beam  density  necessary  to  drive  the  wave  unstable,  n^.  From  Eq.  (1)  a 
marginal  stability  criterion  can  be  established.  For  marginal  stability 


n.of.  +  n  £  +  n,^, 
i  i  e  e  d  b 


where  nj  are  the  densities  of  the  species  and-^j  is  the  power  transferred 
per  particle .  One  then  obtains 


Xe  +  n/±/ne 


critical 
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and 


Y  = 


(4) 


where  Eq.  (4)  is  the  equation  for  the  growth  rate.  If  n^  is  known  then 

the  code  calculates  the  power  transfer  per  particle  for  the  beam,  <^,  and 

the  wave  energy  for  the  wave  being  studied  to  produce  the  growth  rate. 

Obviously  if  n^/ne  <  (n^/n^^^  t^e  mo<^e  damped.  The  power  transfer 
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function  is  calculated  from  the  following  formulas ; 
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For  Maxmillian  Species 
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where  1(b)  is  the  modified  Bessel  function,  b  *  y  M^-j  and  v  =  I  2T  ,m. 

\cj  /  j  \  j  J/ 

As  can  be  seen  from  Eq.  (7)  in  order  to  perform  the  calculations 
one  needs  to  know  the  distribution  function  for  the  beam,  Fb>  One  could 
assume  a  monoenergetic  beam  of  particles  streaming  down  the  field  lines 
but  while  such  a  distribution  will  produce  unstable  waves  it  does  not 
appear  realistic. 

Noting  that  the  debris  streaming  down  the  field  lines  has  already 
interacted  with  the  coupling  shell  through  a  variety  of  instabilities, 
i.e.,  loss  cone  and  mirror  instabilities,  one  needs  a  more  sophisticated 
model.  The  model  must  also  consider  the  effects  of  velocity  dispersion 
because  the  excitation  of  the  scattering  wave  can  occur  at  considerable 
distances  from  the  burt,t  point.  For  a  simple  model  of  these  effects  one 
can  use  the  Fokker-Planck  Equation  to  obtain  a  distribution  which  will 
represent  the  fact  that  the  debris  particles  have  been  scattered  in 
velocity  space.  Inclusion  of  a  source  term  in  the  Fokker-Planck  equation 


1/2 
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allows  one  to  model  the  velocity  dispersion.  That  is,  particles  are 
introduced  into  velocity  space  at  a  given  location.  The  results  of 
running  the  Fokker-Planck  section  of  the  code  is  shown  in  Figure  1. 

g 

The  beam  was  injected  at  200  keV  or  a  velocity  of  1.18  x  10  cm/ sec. 

The  distribution  was  allowed  to  evolve  until  an  appropriate  spread  in 
velocity  was  obtained.  As  can  be  seen  in  Figure  1  the  resulting  distribu¬ 
tion  function  the  particles  have  a  pitch  angle  ratio  of  about  10  to  1, 

Vj  |  to  v^,  consistent  with  the  mirror  ratio  of  the  magnetic  bubble.  In 
addition,  the  particle  density  is  weighted  toward  the  high  velocity  side 
which  is  the  effect  produced  by  velocity  dispersion  along  the  magnetic 
field  line.  The  slow  edge  of  the  distribution  has  a  parallel  velocity 
near  7.3  x  10^  cm/s ec. 

This  distribution  while  not  having  the  exact  details  of  the 
debris  distribution  function  appears  to  have  the  quantitative  details 
which  are  of  interest.  In  the  calculations  to  be  discussed,  this  is  the 
distribution  employed.  It  represents,  we  feel,  a  conservative  estimate 
due  to  the  smoothness  and  the  way  we  introduced  the  beam.  Only  one- 
quarter  of  the  actual  beam  density  was  obtained  so  the  actual  gradients 
are  in  fact  smaller  than  might  be  expected. 

It  should  be  noted  that  the  effect  of  the  velocity  dispersion  is 
very  important  because  the  instability  grows  in  the  positive  slope  region 
of  the  instability.  This  characteristic  has  a  definite  effect  in 
determining  where  the  velocity  will  be  excited  and  what  part  of  the  total 
debris  spectrum  will  be  affected.  These  effects  will  be  discussed  later 


in  this  report. 


Upon  running  the  linear  code  it  was  found  that,  while  the 


fundamental  electrostatic  mode  of  the  ion  cyclotron  wave  could  be 
excited,  very  little  energy  was  transferred  into  this  mode.  Therefore 
growth  rates  were  quite  small  and  insufficient  for  the  anticipated  pulse 
width  of  the  debris.  The  details  of  the  dispersion  relation  for  the  ion 
cyclotron  wave  and  its  higher  harmonics,  the  ion  Bernstein  waves,  can  be 
found  in  Reference  1. 

It  was  found  that  the  ion  Bernstein  waves  could  be  excited  with 
quite  large  growth  rates.  These  rates  approached  the  real  frequency  of 
the  mode  when  n^  was  taken  to  be  20%  of  n^ .  The  parameters  used  for 
these  calculations  are  shown  in  Table  I .  Two  harmonics  of  the  wave  were 
considered.  In  the  first  wave  to  be  studied  the  second  cyclotron  harmonic 
of  the  beam  cyclotron  frequency  resonantes  with  the  wave  near  the  second 
harmonic  of  the  plasma  cyclotron  frequency,  i.e.,:  £u)cb  ~  where 

=  Si  =  2.  In  the  case  of  exciting  the  third  harmonic  of  the  ion 
Bernstein  mode  the  resonate  match  was  =  2,  i  =  3. 

Both  waves  were  found  to  be  easily  excited  because  the  resonant 
nature  of  the  modes  minimizes  the  Landau  and  cyclotron  damping  of  the 
modes.  There  was  considerable  energy  transfer  from  the  beam  to  the  wave 
and  the  growth  rates  approached  the  limits  of  the  code.  Table  II 

show  typical  results  of  the  second  harmonic  calculation.  It  can  be  seen 
from  Table  II  that  the  typical  scale  for  the  wavelengths  are  about  6  km 
in  the  perpendicular  direction  and  30  km  in  the  parallel  direction.  The 
group  velocities  appear  to  be  slow  enough  so  that  the  excitation  process 
will  occur  on  time  scales  faster  than  the  energy  transport.  In  summary 


8 


the  results  from  the  linear  calculations  indicate  that  waves  of  the 
ion  Bernstein  type  will  be  excited  and  so  will  some  of  the  higher 
frequency  waves . 

QUASILINEAR  CALCULATION 

In  this  section  the  quasilinear  estimates  for  the  pitch  angle 
scattering  of  the  streaming  debris  are  addressed.  The  purpose  of  this 
discussion  is  to  ascertain  whether  or  not  quasilinear  effects  would  be 
important.  It  should  be  noted  from  the  outset  that  to  perform  this 
calculation  with  complete  accuracy  one  needs  to  do  the  problem  numerically 
in  two  dimensions.  However,  the  approximation  to  one  dimensional  diffusion 
is  reasonably  good  as  the  perpendicular  electron  field  of  the  wave,  E^, 
is  considerably  larger  than  the  parallel  E|  |  *  Ei  >:>  E|  j  •  The  general 
approach  for  quasilinear  calculations  can  be  found  in  Reference  12. 

As  mentioned,  the  waves  being  considered  are  electrostatic  in 
nature.  This  allows  one  to  write  the  equation  for  the  quasilinear 
behavior  of  the  beam  distribution  as 


where  is  the  Bessel  function  and  ^  is  the  electrostatic  potential. 
For  the  unstable  modes  the  parallel  wave  electric  field  is  sufficiently 
small  that  the  quasilinear  diffusion  reduces  to  a  diffusion  in  v  at  a 
constant  vi i  given  by 
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Equation  (10)  defines  the  quasilinear  time,  i.e.,  the  time  scale  upon 

which  one  expects  quasilinear  affects  to  act.  The  sum  over  the  resonant 

modes  is  related  to  the  sum  over  the  fluctuation  spectrum  through  the 

correlation  time  given  by  Tc  =  l/Akj|V^|  =  (Av^  |  | /v^  |  |  <$co)  where  Av^  |  | /v^ 

is  the  relative  width  of  the  parallel  velocity  distribution  and 

6co  =  w  -  . .  This  time  scale  is  a  measure  of  the  time  the  beam  will 

ci 

stay  in  resonance  with  the  wave.  Although  the  time  T£  can  be  long  com¬ 
pared  to  the  wave  period,  one  finds  that  the  quasilinear  conditions 
remain  valid  for  the  noise  levels  estimated  here  for  many  e-fo]ding  times. 

One  can  define  the  ratio  of  rf  wave  energy  density  to  the  thermal 
plasma  density  by 
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w. 


rf  3 


4  (n  T  +  n,T.\ 
2  \  e  e  x  1/ 


where  the  wave  energy,  W^,  is  defined  in  Reference  1  as 
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where  b 


/2 . 


(12) 


(13) 


For  this  case  E,  =  E  and  it  is  assumed  that  T  =  T.  =  T  for  the  rest  of 
it  x  e  i 

the  calculation.  Substituting  Eq .  (13)  into  Eq .  (12)  returns  the 
following  expression 


Frf  =  Tl  ^  ^ 
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cb 
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(14) 


Using  Eqs.  (10)  and  (11)  as  well  as  the  correlation  time  leads  to  the 
expression  for  the  quasilinear  time 
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From  Eqs.  (14)  and  (15)  one  obtains  the  final  form  for  the  quasilinear 
time 
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where 
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Vb||/Vbi  Avb 


measures  the  localization  of  the  fast  ion  distribution  velocity  space. 

It  is  assumed  that  the  fluctuation  potential  begins  from  thermal  levels 

-3  3 

therefore  (thermal)  =  (2tt)  T  A  k  and  using  Eq.  (13)  the  electric 
field  at  the  thermal  level  is  estimated  to  be 
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One  last  relation  must  be  determined,  that  is  the  growth  time  necessary 

for  the  mode  to  grow  to  nonlinear  or  saturated  levels.  This  can  be 

estimated  by  recalling  that  trapping  effects  will  occur  when  the  noise 

level  has  grown  sufficiently.  Typically  this  would  occur  when  the  trap- 

1/2 

ping  velocity  (e^/m^)  becomes  equal  to  Av^ | |  or  equivalently  when  the 

1/2 

resonant  fast  ion  bounce  frequency,  k||(e<{)/m^)  ,  becomes  equal  to  the 

decorrelation  rate  1/TC*  Applying  this  condition  to  Eq .  (14)  one  arrives 


at  the  relation 


rrf(trap) „  ^  buVL^Y 


Vw: 


with  cg  defined  as  the  sound  speed.  At  this  point  one  can  estimate  the 
time  it  will  take  the  mode  to  encounter  nonlinear  effect  with  the  relation. 


„  J_  o  f  Frf(trap) 
Ts  2y  nlrrf (thermal) 


Using  the  values  contained  in  Tables  I,  II,  and  III,  Eqs.  (16)  and  (19) 
can  be  evaluated.  The  quasilinear  time  estimate  and  the  trapping  time 
estimate  can  be  compared  to  one  another  and  with  the  correlation  time  to 
determine  if  the  quasilinear  effects  have  been  self-consistently  estimated. 

Equations  (16)  and  (19)  along  with  the  relation  for  the  correla¬ 
tion  time  were  numerically  evaluated.  The  correlation  time,  x  ,  was 
found  to  be  in  the  range  of  0.025  to  0.1  sec.  This  implies  that  the 
positive  slope  of  the  distribution  function  such  as  shown  in  Figure  1 
will  be  resonant  with  the  wave  for  times  of  order  0.1  sec.  The  quasi¬ 
linear  time,  x  ,  which  is  the  measure  to  determine  if  strong  pitch  angle 

scattering  will  occur  on  a  time  scale  comparable  with  or  faster  than  the 

_3 

correlation  time,  is  found  to  be  on  order  10  sec  after  the  perturbed 
potential  has  grown  for  'UJ.OS  to  0.1  sec.  The  higher  harmonic  modes 
reach  the  point  of  rapid  quasilinear  diffusion  in  shorter  periods. 

Finally,  the  saturation  time,  x  ,  is  estimated  to  be  approximately  0.05 
to  0.1  sec. 
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The  results  of  evaluating  the  relations  for  x  ,  T  and  T  is 

°  c  q  s 


that  Tc  and  Xfi  operate  on  comparable  time  scales  and  that  on  these  time 


scales  the  wave  potential  has  grown  to  sufficient  levels  to  produce 
quite  rapid  quasilinear  effects  as  reflected  by  the  small  values  of  X 


at  times  on  the  order  of  x  or  x  , 

c  s 


In  Figure  2,  the  results  of  evaluating  the  equation  below  are  shown 


v  (x)  =  v  (0)  + 


E  (eYt  -  1) 
to 


(20) 


where  E^q  is  the  thermal  level  of  the  perpendicular  electric  field.  In 

this  figure  plots  are  made  for  y  ^  2tJci  and  y  'V  3^^ .  They  indicate  that 

v^(t)  will  reach  magnitudes  large  enough  to  remove  the  debris  from  the 

loss  cone  region  on  the  time  scale  of  the  estimated  saturation  time,  X  . 

s 

From  Figure  2,  it  can  be  seen  that  although  the  third  harmonic  wave  begins 
at  a  lower  fluctuation  level,  due  to  its  faster  growth  rate,  it  reaches 
levels  of  sufficient  size  to  cause  the  pitch  angle  scattering  more  rapidly 
than  the  second  harmonic.  It  is  expected  that  waves  with  faster  growth 
rates  would  continue  this  trend. 

Both  the  second  and  third  harmonic  waves  appear  in  Figure  2  to 
produce  sufficient  pitch  angle  scattering  to  explain  the  second  brighten¬ 
ing  phenomena.  The  estimates  indicate  that  the  time  scales  for  the  wave 
to  grow  and  reach  at  least  saturation  amplitudes  are  faster  than  the 


pulse  width  of  debris. 


DISCUSSION 


Due  to  the  very  complicated  relations  between  the  wave  dispersion 
relations  and  the  beam  distribution  function  in  this  calculation  and  the 
lack  of  knowledge  about  the  details  of  the  time  evolution  of  the  debris 
spectrum,  it  is  impossible  to  make  statements  concerning  the  exact  region 
of  the  debris  distribution  function  which  will  excite  the  waves. 

Estimates  can  be  made,  however,  which  will  indicate  probable  areas  of 
interest . 

Before  these  estimates  can  be  made  some  details  of  the  wave 
dispersion  relation  must  be  presented.  It  was  found  during  the  numerical 
calculation  that  the  frequency  at  any  k^  valid  with  the  parameter  regime 
of  the  calculation  was  independent  of  k| | .  It  was  also  found  in  this 
same  parameter  regime  that  the  frequency  was  slowly  varying  in  k^. 
Therefore,  to  good  accuracy  over  the  full  range  of  possible  v^  |  |  one 
could  consider  the  frequency  shift,  u)  -  &ooci ,  to  be  a  constant. 

During  the  course  of  the  calculation  it  was  found  that  the  modes 
that  could  be  excited  had  a  very  limited  range  in  k | |  and  that  this 
preferred  parallel  wavelength  was  about  30  km.  To  some  extent  this 
particular  wavelength  is  a  function  of  the  distribution  used  but  the 
narrowness  of  the  spectrum  is  probably  independent  of  it.  As  can  be 
seen  in  Figure  1,  there  was  a  fair  range  of  V| |  space  in  which  the 
distribution  had  a  positive  slope  but  the  modes  discussed  here  were  the 
only  ones  undergoing  energy  transfer.  Using  the  resonant  condition  for 
beam  wave  energy  transfer,  |  |  =  (u  ~  |  »  anc*  n°ting  that  the 

scale  height  of  the  atmosphere  may  limit  the  length  of  any  given  mode, 
one  can  estimate  an  upper  bound  on  the  parallel  velocity  which  might 


excite  waves.  With  the  scale  height  ^100  km  and  recalling  the 
w  -  %  constant  the  maximum  parallel  beam  velocity  is  estimated  to 

g 

be  about  2.3  x  10  cm/sec.  At  equatorial  latitudes  the  scale  height 
argument  does  not  apply  and  one  may  find  longer  wavelengths  thus  higher 
parallel  velocities.  One  can  make  additional  estimates  by  noting  that 
the  mode  which  appears  to  be  the  most  unstable  in  this  sample  calculation 
has  a  parallel  wavelength  of  about  30  km  and  that  the  correlation  time, 
where  Tc  ^  ^vb  |  |  ^vb  |  |  *  ^as  a  ran8e  °f  about  0.025  sec  to  0.1  sec. 

With  these  parameters  the  range  of  velocities  which  might  excite  the 
waves  is  found  to  be  about  6  x  10  cm/sec  to  1.2  x  10  cm/sec.  However, 
once  the  waves  have  been  excited  debris  of  all  velocities  would  undergo 
strong  pitch  angle  scattering.  Clearly  particles  with  higher  velocities 
than  those  exciting  the  wave  would  also  be  scattered  by  these  waves  but 
the  proportion  of  particles  scattered  at  any  given  energy  would  be  lower. 
This  occurs  because  the  high  energy  particles  emitted  from  the  coupling 
shell  at  the  earliest  times  would  encounter  low  levels  of  wave  turbulence 
while  high  energy  particles  emit  ted  later  would  see  high  fluctuation  levels. 
The  exact  timing  of  such  occurrences  is  difficult  to  ascertain. 

The  time  of  arrival  for  the  particles  with  velocities  in  the 
7  8 

range  6  x  10  to  2.3  x  10  cm/sec  is  the  southern  conjugate  region  (SCR) 
can  be  easily  estimated.  Using  4000  km  as  a  canonical  distance,  one  finds 
that  these  particles  will  arrive  in  the  SCR  in  the  time  frame  of  1.8  to 
7  sec.  However,  since  these  particles  are  expected  to  excite  the  waves 
they  represent  the  sections  of  the  distribution  which  will  be  most 
drastically  isotropized  by  the  waves.  This  would  remove  as  much  as  50% 
of  the  debris  population  in  this  velocity  range  from  the  loss  cone  and 
would  account  for  the  second  brightening  . 
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While  this  discussion  has  been  directed  toward  specific  set  of 
Bernstein  modes,  one  could  expect  a  spectrum  of  both  varying  wavelengths 
and  different  frequency  harmonics  to  be  excited.  The  higher  harmonic 
modes,  l  >  2 ,  would  be  excited  by  the  faster  portions  of  the  streaming 
debris  and  will  be  scattered  by  these  modes.  The  fluctuation  levels  may 
be  expected  to  be  less  than  the  modes  discussed  at  length  here  due  to 

1 

energy  coupling  difficulties  encountered  by  the  higher  harmonic  modes. 

In  any  case  the  range  of  the  spectrum  will  be  greatly  effected  by  the 
exact  details  of  the  debris  distribution  as  a  function  of  time  and  space. 
As  a  final  comment  estimates  made  from  the  SCORPIO  code  runs  indicate 
that  there  would  be  sufficient  beam  density  to  excite  these  waves  in  the 
parameter  region  discussed. 
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CONCLUSIONS 

The  results  of  the  linear  calculations  indicate  that  the  ion 
Bernstein  modes  near  the  second  and  third  harmonic  of  the  plasma 
cyclotron  frequency  will  be  excited  by  the  streaming  debris.  Modes  of 
higher  frequency  may  also  be  excited  but  were  not  considered  in  these 
calculations . 

Quasilinear  estimates  indicate  that  quasilinear  effects  will 
occur  on  time  scales  shorter  than  the  debris  pulse  width  'x-l  sec.  Using 
simple  estimates  for  the  growth  rates  and  the  amplitude  of  the  perpendic¬ 
ular  electric  field  as  a  function  of  time,  it  is  found  that  considerable 
diffusion  of  the  debris  distribution  function  will  occur  on  time  scales 
comparable  with  the  saturation  time,  i.e.,  M).l  sec. 

The  pitch  angle  scattering  predicted  from  the  quasilinear 
estimates  is  more  than  sufficient  to  remove  a  significant  fraction  of 
the  streaming  debris  particles  from  the  loss  cone.  From  the  estimates 
made  of  the  range  of  possible  beam  velocities  that  might  excite  the  wave 
and  undergo  pitch  angle  scattering  it  was  found  that  a  large  percentage 
of  the  total  debris  distribution  was  susceptible  to  these  instabilities. 

The  estimates  appear  to  be  in  agreement  with  the  data  in  the  SCR. 
The  second  brightening  in  this  region  appeared  to  be  nearly  as  bright  as 
the  first  indicating  something  approaching  equal  partition  of  the  energy 
between  the  two  depositions.  From  these  calculations  most  of  the  debris 
distribution  would  be  susceptible  to  the  instabilities  and  one  then  has 
the  possibility  of  up  to  50%  of  the  debris  being  mirrored.  This  of 
course  would  be  the  upper  bound  and  a  more  realistic  number  probably 
involves  factors  of  two.  An  additional  phenomena  was  observed  during 
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the  SCR  deposition.  It  appeared  on  the  film  that  the  altitude  of  deposi¬ 
tion  between  the  first  and  second  brightening  had  changed.  The  region  of 
luminescence  caused  by  the  second  deposition  appeared  at  a  higher  altitude. 
These  results  would  indicate  that  the  energy  of  the  most  energetic 
particles  in  the  debris  being  deposited  was  less  than  at  the  time  of  the 
first  deposition.  This  again  appears  to  be  consistent  with  the  effects 
of  the  instability.  As  the  waves  are  excited,  the  particles  must  lose 
energy.  This  would  have  the  effect  of  decreasing  the  energy  of  the 
portions  of  the  distribution  which  is  exciting  the  waves.  Typically  the 
higher  energy  particles. 

In  summary,  from  the  linear  and  quasilinear  calculations  it  is 
certain  that  the  ion  Bernstein  waves  can  be  excited.  It  appears  that 
they  can  be  excited  to  sufficient  wave  amplitude  to  cause  strong  pitch 
angle  scattering  and  explain  the  second  brightening  seen  in  the  southern 
conjugate  region  during  the  Starfish  event.  It  also  seems  clear  that 
nonlinear  effects  may  make  the  phenomena  discussed  here  even  more  pro¬ 
nounced  if  flattening  of  the  distributions  occurs .  This  flattening  would 
inhibit  the  excitation  of  these  waves  after  the  isotropization .  The 
result  of  this  would  be  that  upon  subsequent  bounces  the  particles  would 
not  be  scattered  out  of  the  loss  cone  by  these  instabilities  and  would 
then  precipitate  for  the  most  part  upon  return  to  the  SCR,  with  particles 
being  returned  to  the  loss  cone  region  by  loss  cone  instabilities.  Work 
on  a  problem  similar  to  this  is  discussed  in  References  13  and  14  where 
the  authors  have  addressed  the  problem  of  electron  trapping  in  the 
radiation  belts. 
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Table  I 

Input  Parameters  for  Ion  Bernstein  Mode 


Electron  Temperature 
Ion  Temperature 
Magnetic  Field 
Electron  Density 
Mass  of  Beam  Particle 
Beam  Density 
Mass  of  Background  Ions 
Initial  Beam  Velocity 


0.07  eV 
0.07  eV 
0.3  G 

1.0  x  10A/c m3 
27  m 

P 

0.2  n 

e 

16  m 

P 

1.18  x  10^  cm/sec  or  200  keV 
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Table  II 


Results  for  the  second  Harmonic  Excitation 


£ 

= 

2 

Growth  Rate  ,y , 

< 

2w  . 

Cl 

ui  -  . 

cx 

'X, 

O.Olu)  .  -+  0.01 

Cl  Cl 

0)  . 
ci 

1.80  x  10^ /sec 

tu  . 
cb 

1.06  x  102/sec 

0) 

% 

3.56  x  102/sec 

k 

1 . E-6 / cm 

i. 

kM 

'X, 

2 .E-6/ cm 

Vthi 

'X, 

9.15  x  10^  cm/sec 

V 

2.97  x  102  cm/sec 

gi 

V  |  | 

% 

1  x  102  cm/sec 

81  1 
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Table  III 


Parameter  for  Quasilinear  Estimate 


% 

v  ^  fto  ,  /k 

bi  cb  i 

iv»ll 

% 

°-2vb|| 

vbl 

'Xy 

°-lvb|| 

4kll 

% 

3  x  10  2/cm 

Ak 

1 

'X/ 

3  x  10  b/cm 

J,(b)e'b 

% 

(b/Z)Vb  <\-  (b/2 

J£(kivi/t0cb) 

'Xy 

a/2 

A3k 

'Xy 

Ak , , Ak2 

1  1  i 

c 

s 

'Xy 

6.47  x  103  cm/ sec 

v,  1.18  x  10b  cm/sec 


Figure  1 


The  debris  distribution  functions  as  used  in  the  linear 
calculation.  The  velocity  of  injection  is  1.18  x  10®  cm/sec. 
The  velocities  on  the  plot  are  in  dimensionless  units.  The 
normalization  velocity  is  5.59  x  106  cm/sec. 


Vj_(0)  =  7.3=<10  cm/sec 

rc  ~  .025-  . 1  sec 
for  both  modes 


is  '  .1  sec  for  1-2 
'  .05  sec  for  a  =  3 


Time  (sec) 

Figure  2 

The  time  evolution  of  the  perpendicular  velocity  of  the 
streaming  debris.  The  initial  v^  is  assumed  0.1  of  v^M  . 
The  parallel  correlation  time,  tc,  and  Ts,  the  estimated 
time  for  parallel  trapping  effects  are  shown  for  each  mode 
In  both  i  =  2  and  £  =  3  case  it  was  assumed  that  Y  ^  &u>ci. 
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naval  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

05CY  ATTN  COOE  552  W.  HOLER 
OICT  ATTN  CODE  0250  C.  BAGGETT 
OICY  ATTN  COOE  81  R.  EASTMAN 


DIRECTOR 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20575 


OICT 

ATTN  COOE 

4700 

OICY 

ATTN  COOE 

4701 

oicy 

ATTN  CODE 

4780 

OICY 

ATTN  COOE 

7500 

OICY 

ATTN  COOE 

7550 

oicy 

ATTN  COOE 

7580 

oicy 

ATTN  COOE 

7551 

OICY 

ATTN  COOE 

7555 

OICY 

ATTN  CODE 

4750 

OICY 

ATTN  CODE 

4127 

TIMOTHY  P.  COFFEY  C25  CYS 
IF  WCLASS,  1  CY  IF  CLASS) 
JACK  0.  BROWN 
BRANCH  HEAD  050  CYS 
IF  UNCLASS,  1  CY  IF  CLASS) 

HQ  COMM  DIR  BRUCE  WALD 
J.  DAVIS 


E.  MCLEAN 
C.  JOHNSON 


commander 

naval  sea  systems  commas 
WASHINGTON,  D.C.  20562 

OICY  ATTN  CAPT  R.  PITKIN 


COMMANDER 

naval  space  surveillance  system 
DAHLGREN,  VA  22448 

OICY  ATTN  CAPT  J.  H.  BURTON 

OFFICER-IN-CHARGE 
naval  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  HD  20910 
OICY  ATTN  CODE  F51 


strategic  ststems  project  office 

DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20576 
OICT  ATTN  NSP-2141 
OICY  ATTN  NSSP-2722  FRED  WIMBERLY 


NAVAL  SPACE  SYSTEM  ACTIVITY 
P.  0,  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CALIF.  90009 

OICY  ATTN  A.  B.  HAZZARC 
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COhfW'GER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAMLGREN,  VA  22448 

01CT  ATTN  CODE  OP-14  R.  SUTLER 


HEADQUARTERS 

electronic  systems  DIVISION/YSEA 

DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  APB,  MA  01731 
01CT  ATTN  TSEA 


COMMAMMNG  OFFICER 
NAVY  SPACE  SYSTEMS  ACTIVITY 
P.O.  BOX  92960 
WORLOWAY  postal  center 
LOS  ANGELES,  CA.  90009 
01CY  ATTN  COOE  52 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  DCKC  MAJ  J.C.  CLARK 


OFFICE  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 

01CY  ATTN  COC*  465 
01CY  ATTN  CODE  461 
01CY  ATTN  COOE  402 
01CY  ATTN  COOE  420 
01CY  ATTN  COOE  421 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND /DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

01CT  ATTN  DC  MR.  LONG 

cot'fwctrc 

AEROSPACE  DEFENSE  COMMAfC/XPO 
OEPARTf€NT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

01CT  ATTN  XPDQQ 
01CT  ATTN  XP 

AIR  FORCE  GEOPHYSICS  LABORATORY 
01751 


COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT -PATTERSON  AFB,  OH  454J3 
01CY  ATTN  NICD  LIBRARY 
01CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  APSC 
GRIFFISS  AFB,  NY  13441 

01CY  ATTN  DOC  LIBRARY /TSLD 
01CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WOR-DWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 


STRATEGIC  AIR  COMWC/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  XPFS  MAJ  B.  STEPHAN 
01CT  ATTN  ADWATE  MAO  BRUCE  BAUER 


OICY 

ATTN  OPR 

HAROLD  GARDNER 

01CT  ATTN 

oicr 

ATTN  OPR-1  JAMES  C.  ULWICK 

OICY  ATTN 

OICY 

ATTN  LKB 

KEMWETM  S.  W.  CHAMPION 

SAMSO/YA 

OICY 

ATTN  OPR 

ALVA  T.  STAIR 

OICY 

ATTN  PHP 

JULES  AARONS 

P.  0.  BOX  92960 

OICY 

ATTN  PHD 

JURGEN  BUCHAU 

WORLDWAY  POSTAL 

OICY 

ATTN  PHD 

JOHN  P.  MULLEN 

LOS  ANGELES,  CA 
OICY  ATTN 

90009 

YAT  CAPT  L. 


BLACK WELDER 


A F  WEAPONS  LABORATORY 
K1RTLAND  AF.B,  NK  87117 


OICY 

ATTN 

SUL 

OICY 

ATTN 

CA 

ARTHUR  H.  GUENTHER 

OICY 

ATTN 

OYC 

CAPT  J.  BARRY 

OICY 

ATTN 

OYC 

J0M4  M.  KANM 

OICY 

ATTN 

OTT 

CAPT  MARK  A.  FRY 

OICY 

ATTN 

OES 

MAJ  GARY  GANONG 

OICY 

ATTN 

OYC 

J.  JAMMI 

AFTAC 

PATRICK  AF8,  Fl  32925 

01CT  ATTN  TF/MAd  WILEY 
01CT  ATTN  7N 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CT  ATTN  AAD  WADE  HJHT 

oicy  attn  aad  allen  johnson 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  t  ACQ 
DEPARTMENT  of  the  air  force 
WASHINGTON,  O.C.  20330 
OICY  ATTN  AFRDQ 

HEADQUARTERS 

electronic  SYSTEMS  D1VISICN/XR 
DEPARTMENT  CP  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  XR  U.  OEAS 


SAMSO/SK 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  AMGELES,  CA  90009 

01CT  ATTN  SKA  (SPACE  COMM  SYSTEMS) 

SAMSO/MN 

NORTON  AFB,  CA  92409 
CMINUTEMAN) 

OICY  ATTN  MNNL  LTC  KEMtCDY 
COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP  A.  LORENTZEN 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.  0.  BOX  5400 
ALBUQUERQUE,  MM  87115 

OICY  ATTN  DOC  CON  FOR  0.  SHERWOOO 

department  of  energy 

LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  OOC  CON  FOR  A.  LABOWITZ 


l.  CLAVIN 
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EGSG,  INC. 

LOS  ALAMOS  DIVISION 

P.  0.  80X  809 

LOS  alamos,  NM  85544 

01CY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  CXX  CON  FOR  TECH  INFO  DEPT 
01CY  ATTN  DOC  CON  FOR  L-589  R.  OTT 
01CY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
01CY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

P.  0.  BOX  1663 

LOS  ALAMOS,  NM  87545 

01CY  ATTN  OOC  CON  FOR  J.  WOLCOTT 
OICY  ATTN  OOC  CON  FOR  R.  F.  TASCHEX 
01CY  ATTN  DOC  CON  FOR  E.  JONES 
OICY  ATTN  DOC  CON  FOR  J.  MALIK 
OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 
OICY  ATTN  DOC  CON  FOR  J.  ZIWV 
OICY  ATTN  DOC  CON  FOR  P.  KEATON 
OICY  ATTN  DOC  CON  FOR  0.  WESTERVELT 


NATIONAL  OCEANIC  i  ATMOSPHERIC  AOMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

OICY  ATTN  R.  GRUBB 

OICY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 
AEROSPACE  CORPORATION 


0.  BOX  92957 
AN&ELES,  CA 

90009 

OICY 

ATTN 

I.  GARFUMCEL 

OICY 

ATTN 

T.  SALMI 

OICY 

ATTN 

V.  JOSEPHSON 

OICY 

ATTN 

S.  BOWER 

OICY 

ATTN 

N.  STOOCWELL 

OICY 

ATTN 

0.  OLSEN 

OICY 

ATTN 

J.  CARTER 

OICY 

ATTN 

F.  MORSE 

OICY 

ATTN 

SMFA  FOR  PWW 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 


SAND  I A  LABORATORIES 
P.  0.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  OOC  CON  FOR  J.  MARTIN 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUGH 

OICY  ATTN  OOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAH.GREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  D1V 


SAM} I A  LABORATORIES 
LIVERMORE  LABORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 
OICY  ATTN  DOC  CON  FOR  T.  COOK 


BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.  0.  BOX  983 
BERKELEY,  CA  94701 

OICY  ATTN  J.  WORKMAN 

BOEING  COMPANY,  THE 
P.  0.  BOX  3707 
SEATTLE,  WA  98124 

OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KEWEY 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
IPAPS,  B-019 
LA  JOLLA,  CA  92093 

OICY  ATTN  HENRY  G.  BOOKER 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  OOC  CON  FOR  0.  GALE 


BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 

OICY  ATTN  ROMEO  A.  DELIBERIS 


OTHER  GOVERNMENT 


CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RD/SI,  RM  5G48,  HQ  BLOG 
WASHINGTON,  D.C.  20505 

OICY  ATTN  OSI/PSID  RM  5F  19 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANOARDS 
WASHINGTON,  D.C.  20234 

CALL  CORRES:  ATTN  SEC  OFFICER  FOR) 
OICY  ATTN  R.  MOORE 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  8  INEO  ADMIN 
BOULDER,  CO  80303 

OICY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
OICY  ATTN  W.  UTLAUT 
OICY  ATTN  0.  CROMBIE 
OICY  ATTN  L.  BERRY 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

OICY  ATTN  D.  B.  COX 
OICY  ATTN  J.  P.  GILMORE 

COMPUTER  SCIENCES  CORPORATION 
6565  ARLINGTON  BLVD 
FALLS  CHURCH,  VA  22046 
OICY  ATTN  H.  BLANK 
OICY  ATTN  JOHN  SPOOR 
OICY  ATTN  C.  NAIL 


COMSAT  LABORATORIES 
LI NTH  I CUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.  T.  FARLEY  JR 
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I-A'WW'P  ’***■*!*" 


ELECTROSPACE  SYSTEMS,  INC. 

BOX  1 359 

RICHARDSON,  TX  75080 

01CY  ATTN  H.  LOGSTON 

01CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL  INC. 

495  JAVA  DRIVE 

SUNNYVALE,  CA  94086 

01CY  ATTN  J.  ROBERTS 
01CY  ATTN  JAF€S  MARSHALL 
01CY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  t  COMMUNICATIONS  CORP 

5939  FABIAN  WAY 

PALO  ALTO,  CA  94503 

OlCY  ATTN  J.  T.  MATTINGLEY 


GENERAL  ELECTRIC  C OFF* ANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GOOOARD  8LVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OlCY  ATTN  M.  H.  80RTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

OlCY  ATTN  F.  REI8ERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUOIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OlCY  ATTN  DAS I AC 
OlCY  ATTN  DON  CHAM3LER 
OlCY  ATTN  TOM  BARRETT 
OlCY  ATTN  TIM  STEPHANS 
OlCY  ATTN  WARREN  S.  KNAPP 
OlCY  ATTN  WILLIAM  MCNAMARA 
OlCY  ATTN  B.  GAMBILL 
OlCY  ATTN  MACK  STANTON 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

OlCY  ATTN  G.  MILLMAN 


GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  93111 

OlCY  ATTN  JOHN  ISE  JR 
OlCY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
OlCY  ATTN  T.  N.  DAVIS  (UNCL  0N.Y) 

OlCY  ATTN  NEAL  BROWN  (UNCL  0N.Y) 

OlCY  ATTN  TECHNICAL  LIBRARY 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NCEDHAM,  MA  02194 

OlCY  ATTN  MARSHAL  CROSS 


ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
UR6ANA,  1L  61803 

OlCY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
URBANA,  IL  60680 

(ALL  COR  RES  ATTN  SECURITY  SUPERVISOR  FOR) 
OlCY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 

OlCY  ATTN  J.  M.  AE1N 
OlCY  ATTN  ERNEST  BAUER 
OlCY  ATTN  HANS  WOLF  HARD 
OlCY  ATTN  JOEL  BENGSTON 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OlCY  ATTN  DONALD  HANSEN 

INTL  TEL  C  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OlCY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  C AMI  NO  DEL  MAR 
DEL  MAR,  CA  92014 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  M3  20810 

OlCY  ATTN  DOCUMENT  LIBRARIAN 
OlCY  ATTN  THOMAS  POT  EMI  A 
OlCY  ATTN  JOHN  DASSOULAS 

LOCKHEED  MISSILES  C  SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  94088 

OlCY  ATTN  DEPT  60-12 
OlCY  ATTN  D.  R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OlCY  ATTN  MARTIN  WALT  DEPT  52-10 
OlCY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OlCY  ATTN  W.  L.  IFHOF  DEPT  52-12 

KAMAN  SCIENCES  CORP 
P.  0.  BOX  7463 
COLORADO  SPRINGS,  CO  80933 
OlCY  ATTN  T.  MEAGHER 

lincabit  corp 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OlCY  ATTN  IRWIN  JACOBS 

M.I.T.  LINCOLN  LABORATORY 
P.  0.  BOX  75 
LEXINGTON,  MA  02173 

OlCY  ATTN  DAVID  M.  TOWLE 
OlCY  ATTN  P.  WALDRON 
OlCY  ATTN  L.  LOUGHLIN 
OlCY  ATTN  D.  CLARK 

martin  marietta  corp 

ORLANDO  DIVISION 
P.  0.  BOX  5837 
ORLANDO,  Fl  32805 

OlCY  ATTN  R.  HEFFNER 
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MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  8EACH,  CA  92647 
01CT  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MRO Z 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.  W,  HALPR1N 
01CY  ATTN  TEOMCAL  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  9J101 
01CY  ATTN  P.  FISCHER 
01CY  ATTN  «.  F,  CREVIER 
01CY  ATTN  STEVEN  L.  GUTSCHE 

oicr  attn  o.  Sapped i elo 

01CY  ATTN  R.  BOGUSCH 

01CY  ATTN  R.  HEM3RJCK 

01CY  ATTN  RALPH  KILS 

01CY  ATTN  OAVE  SOWLE 

01 CY  ATTN  F .  FA  JEN 

01CY  ATTN  M.  SCHE1BE 

01CY  ATTN  CONRAD  L.  LONGMIRE 

01CY  ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  TN£ 

P.  0.  BOX  208 
BEDFORD,  MA  01750 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.  E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BCVO 
MCLEAN,  VA  22101 

01CY  ATTN  W.  HALL 
01CY  ATTN  w.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERF I ELD  BLVD. 

SANTA  MONICA,  CA  901*09 

01CY  ATTN  E.  C.  FIELD  JR 

FENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
518  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

FHOTOFCTR1CS,  INC. 

942  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.  J.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  94701 

OICY  ATTN  A.  THOMPSON 


R  t  0  ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WRIGHT  JR 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 


RAF®  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZ I  AN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  Et®  AVEIWJE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2351 
LA  JOLLA,  CA  92038 

01CT  ATTN  LEWIS  M.  L1NSON 
01CT  ATTN  DANIEL  A.  HAMLIN 
01 CT  ATTN  D.  SACHS 
OICY  ATTN  E.  A.  STRAKER 
01CT  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCOOUGALL 

RAYTHEON  CO. 

528  boston  post  road 
SUDBURT,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

science  applications,  inc. 

HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUVTSVILLE,  Al  55805 

OICY  ATTN  DALE  H.  DIVIS 

SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WESTPARX  DRIVE 
MCLEAN,  VA  22101 

OICY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  DRIVE 
PLEASANTON,  CA  94566 
OICY  ATTN  SZ 


SRI  INTERNATIONAL 
533  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

OICY  ATTN  DONALD  NEILSON 
01CT  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  l.  L.  COBB 
OICY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  R1NO 
01CT  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  ATTN  RAY  L.  LEADABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICT  Attn  j,  PETERSON 
01CT  ATtn  R.  HAKE,  JR. 

OICT  ATTN  V.  GONZALES 
OICY  attn  D.  MCDANIEL 
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